Photon storage in A-type optically dense atomic media. IV. Optimal control using 

gradient ascent 
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We use the numerical gradient ascent method from optimal control theory to extend efficient 
photon storage in A-type media to previously inaccessible regimes and to provide simple intuitive 
explanations for our optimization techniques. In particular, by using gradient ascent to shape 
classical control pulses used to mediate photon storage, we open up the possibility of high efficiency 
photon storage in the non-adiabatic limit, in which analytical solutions to the equations of motion do 
not exist. This control shaping technique enables an order-of-magnitude increase in the bandwidth 
of the memory. We also demonstrate that the often discussed connection between time reversal 
and optimality in photon storage follows naturally from gradient ascent. Finally, we discuss the 
optimization of controlled reversible inhomogeneous broadening. 
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I. INTRODUCTION 

Faithful mapping between quantum states of light (fly- 
ing qubits) and quantum states of matter (storage and/or 
memory qubits) is an important outstanding goal in the 
field of quantum information processing and is being pur- 
sued both theoretically and experimentally by a large 
number of research groups around the world. Photon 
storage in A-type atomic media is a promising avenue for 
achieving this goal. In a recent theoretical paper [l|, we 
unified a wide range of protocols for photon storage in A- 
type media, including the techniques based on Electro- 
magnetically Induced Transparency (EIT), off- resonant 
Raman interactions, and photon-echo. In Ref. we 
also demonstrated equivalence between all these proto- 
cols and suggested several efficiency optimization proce- 
dures, some of which have since been demonstrated ex- 
perimentally J31. In the three preceding papers of this 
series, Refs. [IJl, which we will refer to henceforth as 
papers I, II, and III, we presented some details and many 
extensions of the analysis of Ref. [l| ■ Most of the results 
in Ref. fl'] and in papers I, II, and III were obtained based 
on physical arguments and on exact solutions available in 
certain limits. However, the optimization problems dis- 
cussed there fall naturally into the framework of optimal 
control problems, for which powerful numerical optimiza- 
tion methods exist ^,'7*1 . Thus, in the present paper, we 
apply these optimal control methods to the problem of 
photon storage. As a result, we open up the possibil- 
ity of efficient photon storage in previously inaccessible 
regimes by increasing the bandwidth of the memory and 
provide simple intuitive understanding for the optimiza- 
tion methods underlying photon storage. 

We refer the reader to paper I for a comprehensive 
introduction to photon storage in A-type atomic media 
and for the full list of references. Here we summarize 



only a few important points. In a typical photon storage 
protocol, an atomic ensemble with A-type level structure 
shown in Fig. [1] is assumed to start with all N atoms 
pumped into the metastable state \g). The incoming 
quantum light mode is coupled to the \g) — \e) transition 
with a collectively enhanced coupling constant g\/N and 
is mapped onto the collective coherence (called a spin 
wave) between the metastable states \s) and \g) using 
a classical two-photon-resonant control pulse with time- 
dependent Rabi frequency U,(t). Ideal mapping of the 
light mode onto the spin wave and back can be achieved 
in an ensemble that has infinite resonant optical depth d 
on the \g) — |e) transition. However, despite the existence 
of proposals for achieving high values of d Q , in most cur- 
rent experiments d (or the cooperativity parameter C for 
ensembles enclosed in a cavity [9]) is limited to d ~ 10 due 
to experimental imperfections such as competing four- 
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FIG. 1: (Color online) A-type medium coupled to a quantum 
field (dashed) with a collectively enhanced coupling constant 
qVN and a two-photon-resonant classical field (solid) with 
time-dependent Rabi frequency n(t) . 
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wave mixing processes [l^, spatially- varying light shifts 
0, number of atoms in a trap [111, [l^], or inhomoge- 
neous broadening and short interaction lengths [13, [Tj] . 
As a result of the limited optical depth, the experimen- 
tally demonstrated efficiencies for the light-matter inter- 
face are low, which makes the optimization of photon 
storage protocols at finite values of d crucial. The opti- 
mization in Ref. in papers I, II, and III, as well as in 
the present paper relies on the knowledge of the shape 
of the incoming photon mode. Note that such knowl- 
edge is not incompatible with storing unknown quantum 
states because the mode usually acts simply as a carrier 
while the information is stored in the quantum state of 
the harmonic oscillator corresponding to this mode [l5l |. 
A different type of problem is the storage of an unknown 
mode or, equivalently, the storage of multiple photonic 
modes within an ensemble [16j. While we believe that 
the optimization procedures considered here will proba- 
bly also be relevant to this situation, we shall not discuss 
it in more detail here. 

The main tool used in this paper is a numerical it- 
erative optimization with respect to some set of control 
parameters, which are updated to yield higher photon 
storage efficiency at every iteration. Such iterative opti- 
mization methods are a standard tool in applied optimal 
control theory 0| • These methods and their variations 
are already being used in a variety of applications includ- 
ing laser control of chemical reactions [H, , design 
of NMR pulse sequences [1^, loading of Bose-Einstein 
condensates into an optical lattice [2l|, atom transport 
in time-dependent superlattices [23] , quantum control of 
the hyperfine s pin of an atom [23] , and design of quan- 
tum gates [l^, llHl • Although advanced mechanisms for 
updating the control parameters from one iteration to 
the next exist and exhibit superior convergence charac- 
teristics [6, 21, 26, 27, 28], we will concentrate in the 
present paper on optimization via a simple gradient as- 
cent method B S 0, [23], except for Sec. UHl where ad- 
vanced updating mechanisms will also be used. Gradi- 
ent ascent methods are often more efficient than simple 
variations of the control parameters using, e.g., genetic 
algorithms. Moreover, we will show that gradient ascent 
optimization has the advantage that it can often be un- 
derstood physically and can provide deeper intuition for 
the photon storage problem. In particular, in papers I, II, 
and III, we used involved physical arguments and exact 
analytical solutions available in certain limits to derive a 
time-reversal-based iterative optimization with respect to 
the shape of the incoming photon mode. In the present 
paper, we show that these time-reversal iterations and 
the general and often discussed connection between opti- 
mality and time reversal in photon storage p]. [29l[30ll33j 
naturally follow from the gradient ascent method. The 
results of papers I, II, and III are, however, still cru- 
cial since they show in certain cases that the solutions 
obtained via the local gradient ascent method represent 
global, rather than local, optima. 

In addition to considering optimization with respect to 



the shape of the input mode, we consider in the present 
paper optimization with respect to the storage control 
field. In particular, we show that shaping the control 
field via the gradient ascent method allows for efficient 
storage of pulses that are an order of magnitude shorter 
than when the control field is optimized in the adiabatic 
approximation discussed in Ref. [1] and in papers I and 
II. In other words, this new control shaping method in- 
creases the bandwidth of the memory. Finally, we discuss 
the performance of optimal control pulses in the con- 
text of photon storage via controlled reversible inhomo- 
geneous broadening (CRIB) [s^]- In particular, assum- 
ing one is interested in storing a single known incoming 
light mode and assuming one can shape control pulses 
with sufficient precision, we are not able to identify any 
advantages of CRIB-based photon storage compared to 
photon storage with optimal control pulses in homoge- 
neously broadened media. 

The remainder of the paper is organized as follows. 
In Sees, ini mil and IIVI we show how gradient ascent 
can be used to optimize with respect to the control field, 
the input mode, and the inhomogeneous profile, respec- 
tively. We summarize the discussion in Sec. |V| and, fi- 
nally, present some details omitted in the main text in 
the Appendixes IXllDl 



II. OPTIMIZATION WITH RESPECT TO THE 
STORAGE CONTROL FIELD 

In principle, both the incoming light mode and the 
classical control pulse may be adjusted to maximize the 
light storage efficiency. However, it is often easier to 
vary the classical control pulse. In particular, the pho- 
tonic state we wish to store may be some non-classical 
state generated by an experimental setup, where we can- 
not completely control the shape of the outgoing wave 
packet. This is, e.g., the case for single photons gener- 
ated by parametric down conversion |16l . 1311 or by sin- 
gle nitrogen- vacancy centers in diamond [32], where the 
shape of the wave packet will be, respectively, set by the 
bandwidth of the setup and the exponential decay associ- 
ated with spontaneous emission. Alternatively, the wave 
packet may also be distorted in an uncontrollable way by 
the quantum channel used for transmitting the photonic 
state [isj . In this section, we therefore discuss optimiza- 
tion with respect to the storage control field in both the 
cavity model (Sec. Ill A|) and the free space model (Sec. 

MM- 



A. Cavity model 

As discussed in papers I and II, the cavity model, in 
which the atomic ensemble is enclosed in a cavity, is the- 
oretically simpler than the free space model because only 
one collective atomic mode can be excited. In addition, 
as shown in papers I and II, the cavity setup can yield 
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higher efhciencies in certain cases than the free space 
model due to the enhancement of the optical depth by the 
cavity finesse and due to (for certain spin wave modes) 
better scaling of the error with the optical depth d {1/d 
in the cavity vs. 1 /-v/d in free space) . We, therefore, start 
with the cavity model. As in paper I, to get the closest 
analogy to the free-space regime, we will discuss in the 
present paper only the so-called "bad cavity" limit, in 
which the cavity mode can be adiabatically eliminated. 
However, the method of gradient ascent can easily be 
applied outside of this limit, as well. 

To simplify the discussion, we first consider the sim- 
plest example, in which one stores a given resonant input 
mode into a homogeneously broadened ensemble enclosed 
in a cavity and having negligible spin-wave decay rate. 
It is important to note that, because only one spin- wave 
mode is accessible in the cavity model, the retrieval ef- 
ficiency is independent of how the storage is done (see 
paper I). This makes it meaningful to optimize storage 
separately from retrieval [the latter does not have to be 
optimized since its efficiency depends only on the coop- 
erativity parameter (see paper I)] . 

We follow the derivation of paper I to adiabatically 
eliminate the cavity mode and to reduce the equations 
of motion to the following complex number equations on 
the time interval t G [0,T]: 



P{t) = -j{l + C)P{t) 

s{t) = m{t)p{t). 



(2) 



Here the optical polarization P{t) on the \g) — |e) transi- 
tion and the spin polarization S{t) on the \g) — |s) tran- 
sition satisfy initial conditions P{0) — and S{0) = 0, 
respectively, corresponding to the absence of atomic ex- 
citations at t = 0. In this example, the shape of the 
incoming mode £in{t) is assumed to be specified, real, 
and normalized according to dt£^^{t) = 1. 7 is the 
decay rate of the optical polarization and C is the col- 
lectively enhanced cooperativity parameter equal to the 
optical depth of the atomic ensemble times the cavity 
finesse. The goal is to find the slowly varying control 
field Rabi frequency envelope fl{t) (assumed to be real) 
that maximizes the storage efficiency 77s = |'5'(7')p. [To 
avoid carrying around extra factors of 2, fl{t) is defined 
as half of what is usually called the Rabi frequency: it 
takes time 7r/(2r2) to do a tt pulse]. For the moment, 
we suppose that there is no constraint on the energy of 
the control pulse and return to the possibility of includ- 
ing such a constraint below. It is worth noting that due 
to their linearity, the equations of motion (and all the 
results of the present paper) apply equally well both to 
classical input fields with pulse shapes proportional to 
£in(t) and to quantum fields whose excitations are con- 
fined to the mode described by £in{t)- The efficiency ij 
is thus the only parameter required to fully characterize 
the memory (see paper I). 

Since the optimization of ?7s is constrained by the equa- 
tions of motion ([T]) and we introduce Lagrange mul- 



tipliers P{t) and S{t) to ensure that the equations of 
motion are fulfilled, and turn the problem into an uncon- 
strained maximization of 01 

J = S{T)S*{T) 

i-T 

dt 



dt 



P*{ -P-jil + C)P + inS 

in I + C.C. 

S* i-S + inp] +C.C. 



(3) 



where c.c. stands for the complex conjugate [34| . 

The optimum requires that J is stationary with respect 
to any variation in P, S, and Q. As shown in Appendix 
|A| setting J to be stationary with respect to variations 
in P and S requires that the Lagrange multipliers (also 
referred to as the adjoint variables) P and S satisfy the 
equations of motion 



p = -f{i + c)p + ins, 

§ = iQP, 

subject to boundary conditions at time t — T 



P{T) 
S{T) 



0, 

Sir). 



(4) 
(5) 



(6) 
(7) 



These are the same equations as for S and P [Eqs. ([T|) 
and ^] except that there is no input field and that the 
decay with rate 7(1 -I- C) is replaced with growth, which 
will function as decay for backward evolution. This back- 
ward evolution, in fact, corresponds to retrieval with the 
time-reversed control field and can be implemented ex- 
perimentally as such (see papers I, II and Ref. 0]). It 
is satisfying to have obtained this purely mathematical 
and simple derivation of the often discussed connection 
between optimality and time reversal in photon storage 
[E i, 1 i, S, m S m, HI. As explained in the intro- 
duction to paper I and as shown in detail in paper II, this 
connection goes beyond the perfect reversibility of uni- 
tary evolution discussed in Refs. [2^,[33,|3l| by including 
systems with non-reversible dynamics, as exemplified, for 
example, by the decay rate 7 in the present model. 

Eqs. (HI)-® ensure that J is stationary with respect to 
variations in P and S. To find the optimum it remains 
to set to zero the functional derivative of J with respect 
to ri. This functional derivative is given by 



6J 

6n{t) 



-2Im [S*P~PS*] , 



(8) 



where "Im" denotes the imaginary part. 

In general, if one has a real function of several vari- 
ables, one way to find a local maximum is to pick a 
random point, compute partial derivatives at that point, 
move a small step up the gradient, and then iterate. The 
same procedure can be applied to our optimal control 
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problem fd, ^ . The gradient ascent procedure for finding 
the optimal storage control pulse ^l{t) is to take a trial 
il{t) and then iteratively update r2(t) by moving up the 
gradient in Eq. ^ according to 

n{t)-^n{t)~\lin[S*P-PS*]. (9) 
A 

where 1/A regulates the step size. In order to compute 
the right hand side of Eq. one has to evolve the 
system forward in time from t — to t = T using Eqs. 
dll) and (HI) to obtain S{t) and P{t). Then project the 
final atomic state described by S{T) and P{T) onto S 
according to Eqs. (©and Q to obtain P(T) and S{T). 
Then evolve S and P backwards in time from t = T to 
t — according to Eqs. ([4]) and ([5]). 

In general, as in any gradient ascent method, the step 
size 1/A in Eq. ^ has to be chosen not too big (one 
should not go up the gradient so quickly as to miss the 
peak) but not too small (in order to approach the peak 
relatively quickly). To achieve faster convergence, one 
could use a different step size 1/A for each iteration; but 
for the problems considered in the present paper, conver- 
gence is usually sufficiently fast that we do not need to do 
this (unless the initial guess is too far from the optimum, 
in which case changing A a few times helps). Moreover, 
in some optimization problems [l^j, 1/A has to be chosen 
such that it depends on the argument of the function we 
are trying to optimize, i.e., in this case the time i; this 
is not required for the present problems, and 1/A is just 
taken to be a constant. 

For example, let us take C = 1, T7 — 10, and a 
Gaussian-like input mode 

fi,(t) = ^(e-30(t/T-o.5)^ _ e-7-5)/Vr, (10) 

where A « 2.09 is a normalization constant and where 
the mode is chosen to vanish at t — and t = T for com- 
putational convenience. Starting with an initial guess 
n{t) = \/tv^ and using A = 0.5, it takes about 45 it- 
erations for the efficiency to converge to within 0.001 of 
the optimal efficiency of C/(l + C) = 0.5 (see paper I 
for the derivation of this formula). If, however, A is too 
small (e.g. A = 0.1), then the step size is too large, and, 
instead of increasing with each iteration, the efficiency 
wildly varies and does not converge. 

We now compare the optimal control field shaping to 
the adiabatic control field shaping presented in paper I. 
We first take C = 10 and consider the input mode in 
Eq. (fTO|) with T — 50 /j. Following paper I, we calcu- 
late the storage control field using the adiabatic equa- 
tions [Eq. (26) in paper I], then numerically compute the 
storage efficiency with this control field, and multiply it 
by the complete retrieval efficiency C/{1 + C) to obtain 
the total efRciency. Since we are in the adiabatic limit 
(TC7 = 500 » 1), the resulting total efficiency is equal 
to the maximum possible efRciency C^/(l + C)^ = 0.83 
(see paper I). Fig.[2)Ja) shows the input mode in Eq. ITUl) 
(dashed line) and the adiabatic storage control field (dot- 
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FIG. 2: Adiabatic (dotted) and optimal (solid) control fields 
for the storage of a Gaussian- like input mode £in{t) (dashed) 
in the cavity model with C = 10 and T = 5O/7 (a) and 
T = 0.5/7 (b). The four different optimal control pulses 
correspond to four different initial guesses for the gradient 
ascent optimization. The adiabatic control field agrees with 
the optimal one in the adiabatic limit {TC'y 2> 1) (a) and 
deviates from it otherwise (b). 



ted line). The optimal control field shaping using gradi- 
ent ascent via Eq. ^ also yields the maximum possible 
efSciency C^/(l -I- C)^ = 0.83 independent of the initial 
guess for f2(i). The four solid lines in Fig. [D^a) show 
f2(t) resulting from optimal control field shaping for four 
different initial guesses, il{t)/'j = 0.2, 1, 2, and 3. The 
four optimal control fields and the adiabatic control field 
agree except at small times. The reason for the disagree- 
ment is that the dependence of storage efficiency on the 
front section of the control field is very weak because this 
section affects only the front part of the excitation, and 
a large part of this anyway leaks out at the back end of 
the atomic ensemble. In fact, the dependence is so weak 
that gradient ascent leaves the front part of the initial 
guesses almost unperturbed. 

It is worth noting that, in general, gradient ascent 
methods are not guaranteed to yield the global optimum, 
and the iterations may get trapped in a local maximum. 
However, for our photon storage problem, we know what 
the global optimum is in some cases. In particular, we 



have shown in paper I (for the cavity model) and in paper 
11 (for the free space model) that, in the adiabatic limit, 
adiabatic control field shaping yields the global optimum. 
Since control shaping via gradient ascent agrees with the 
adiabatic shaping in this limit, we have a strong indica- 
tion that gradient ascent always yields the global opti- 
mum also outside of the adiabatic limit. The global op- 
timum is here the (unique) maximum possible efhciency, 
which, within the numerical error, is achievable for a va- 
riety of control fields due to the lack of sensitivity to the 
control field for small times (see Fig. [2]). 

We now repeat the same steps except that we use 
T = 0.5/7. The resulting control fields are shown in Fig. 
HJb). Again the four optimal control fields correspond 
to diS^erent initial guesses [^l{t)/^ = 2, 5, 8, and 11]. 
The adiabatic control field now differs from the optimal 
one on the entire time interval. The reason is that the 
adiabatic limit {TCj ^ 1) is not satisfied to a sufficient 
degree {TCj = 5), and, as a result, the adiabatic approx- 
imation does not work well. Indeed, the efficiency yielded 
by the adiabatic control (0.49) is much smaller than that 
yielded by the optimal control (0.81). Physically, the 
breakdown of the adiabatic approximation means that 
the optical polarization P{t) no longer follows the spin 
wave S{t) adiabatically, but rather evolves dynamically 
according to the full differential equation ([T]). Since in 
this regime (TC7 1) the optimal control field is turned 
on abruptly following a time period when it is off [see Fig. 
I^b)], the optimal storage procedure acquires some char- 
acteristics of photon-echo type fast storage [l|, |3, 0, HI] ■ 
In fast storage, the input pulse is first absorbed on the 
|e) — \g) transition in the absence of the control field, and 
is then mapped to the |s) — \g) coherence via a control tt 
pulse. This connection is not surprising since fast stor- 
age is indeed optimal for certain input modes of duration 
T ~ 1/(C7). Finally, we note that all the initial guesses 
for fl that we tried yielded the same optimal control (up 
to the unimportant front part) and the same efficiency, 
which is a signature of the robustness of the optimiza- 
tion protocol and is another strong indication that, for 
this optimal control problem, gradient ascent yields the 
global, rather than local, optimum. 

Having performed the comparison of the control fields 
generated by adiabatic shaping and by gradient ascent, 
we turn to the investigation of the dependence on C and 
on TC'j of the efficiency achieved by these two meth- 
ods. In Fig. [3ja), we compare the efficiency of storage 
followed by retrieval of the input mode of Eq. ([TU)) ob- 
tained using the adiabatic control field (dotted lines) and 
using the control found via gradient ascent (solid lines). 
The efficiencies are plotted as a function of TCj for three 
indicated values of C {— 1, 10, 100). Dashed lines corre- 
spond to C^/(H-C)^, the maximum efficiency possible at 
any given C. We note that the dotted lines have already 
been shown in Fig. 2(a) of paper I. According to the ar- 
guments presented in papers 11 and III, we note that it 
is impossible to retrieve into a mode much shorter than 
1/(7(7), and hence, by time-reversal, it is impossible to 
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FIG. 3: (a) The total efficiency of storage followed by retrieval 
for the Gaussian-like input mode in Eq. pop using adiabatic 
equations (dotted) and gradient ascent (solid) to shape the 
storage control field. Results are shown as a function of TC'j 
for the indicated values of C (=1, 10, 100). The dashed lines 
are CV(1+C)^ the maximum efficiency possible at any given 
C. (b) Same for Sinit) = 1/VT. 



efficiently store such a short mode. Fig. [3)Ja) confirms 
that indeed, when TC7 ^ 1, even optimal controls can- 
not give high efficiency. Using gradient ascent instead of 
adiabatic shaping, one can, however, efficiently store in- 
put modes that are about an order of magnitude shorter 
and, thus, an order of magnitude larger in bandwidth. 
It is worth repeating that although the method of gradi- 
ent ascent is generally not guaranteed to yield the global 
maximum, the fact that it does give the known global 
maximum in the limit TC7 ^ 1 suggests that it proba- 
bly yields the global maximum at all values of TCj. 

To confirm the robustness and generality of the opti- 
mization procedure, we show in Fig. E^b) the results of 
the same optimization as in Fig. [3]Ja) but for a square 
input mode Emit) — 1/VT instead of the Gaussian- like 
input mode of Eq. (fTO]) . As in Fig.[3lja), we see that gra- 
dient ascent control shaping improves the threshold in 
the value of TCj, where efficiency abruptly drops, by an 
order of magnitude. This can again be interpreted as an 
effective increase in the bandwidth of the memory by an 
order of magnitude. The optimal storage efficiency for 
the square input pulse falls to half of the maximum at 
smaller TC7 than for the Gaussian- like input pulse be- 
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cause the latter has a duration (half-width at half maxi- 
mum, for example) significantly shorter than T [see Eq. 
(fTO| or Fig. [2]. On the other hand, as TC7 is increased, 
the maximum is approached slower for the square input 
mode than for the Gaussian-like mode. This is because 
the high frequency components contributed by the sharp 
edges of the square pulse are difficult to store. 

Most experiments have features that go beyond the 
simple model we have just described. Therefore, in Ap- 
pendix [Bj we generalize this model and the optimization 
procedure to include the possibility of complex control 
field envelopes f2(i) and input mode envelopes £in(i), 
nonzero single-photon detuning A and spin wave de- 
cay rate 7s, and (possibly reversible [sO]) inhomogeneous 
broadening. Our model of inhomogeneous broadening is 
applicable both to Doppler broadening in gases and to the 
broadening of optical transitions in solid state impurities 
caused by the differences in the environments of the im- 
purities [3Q|. For the case of Doppler broadened gases, we 
also allow for the possibility of modeling velocity chang- 
ing collisions with rate 7c. Finally, in Appendix |Bl we 
also show how to take into account the possibility that 
the classical driving fields available in the laboratory are 
not sufficiently strong to realize the optimal control fields, 
which may be the case for short input modes and/or large 
single-photon detuning A, both of which require control 
pulses with large intensities. 

Although a comprehensive study of optimization for 
A 7^ is beyond the scope of the present paper, we will 
now prove that the maximum efficiency for A 7^ is 
exactly equal to the maximum efficiency for A = 0. Sup- 
pose we know the control field Oo(i) that achieves the 
optimum for a given resonant input £in{t). Then, for an 
input at A 7^ with the same envelope fin(i), we can 
construct the control field Vl{t) as a sum of two parts 
[written in the two-photon-resonant rotating frame as in 
Eqs. dHIl) and (|B2l) ] 



iAot 



iAt 



(11) 



The first part is a far-detuned control (A2 ^ ^^2,7) that 
Stark-shifts level |e) into resonance with the input (i.e. 
such that f72/A2 = A), while the second part is resonant 
with the Stark-shifted |e) — |s) transition and has the en- 
velope equal to the optimal resonant control. The reason 
why an extra detuning A is needed to bring the second 
term in two-photon resonance is because fl2 Stark-shifts 
both |e) and |s) by A. The resulting efficiency must be 
equal to the optimal resonant efficiency up to an error on 
the order of the small population mixing between |e) and 
|s) caused by Q2', that is, ~ {^2/^2)^ — A/A2. To verify 
mathematically that the control in Eq. (fTTj) works, one 
can write P and as a sum of a slowly varying piece and 
a rapidly oscillating piece, extract separate equations for 
the rapidly and slowly oscillating variables, and finally 
adiabatically eliminate the rapidly oscillating variables. 
We have also numerically verified the performance of the 
control in Eq. (fTTj) and the scaling of the error (~ A/A2) 
by integrating the equations of motion for the case of ho- 



mogeneous broadening at several different values of T7 
and C for the pulse shape in Eq. pU)) . Thus, the op- 
timal off-resonant efficiency is greater than or equal to 
the optimal resonant efficiency for the same input enve- 
lope £in(t). Carrying out the same argument backwards 
[i.e. using il2(0 to shift |e) out of resonance], we con- 
clude that the optimal efficiency must be the same on 
and off resonance. When applying this idea in practice, 
one should, of course, realize that, in addition to a pos- 
sible technical limit on the available control power, the 
three-level approximation and the rotating-wave approx- 
imation may start to break down for sufficiently large 
values of A2 and ^2- 



B. Free space model 



Although the cavity model discussed in Sec. Ill Al ls the- 
oretically simpler and results, in certain cases, in higher 
efficiencies than the free space model, the latter is easier 
to set up experimentally. Moreover, because of the ac- 
cessibility of a large number of spin wave modes, the free 
space model can provide higher efficiencies in some other 
cases (see paper II) and can, in principle, function, as a 
multi-mode memory. Therefore, we turn in the present 
section to the analysis of the free space model. 

To demonstrate how optimization with respect to the 
control field works in the free space model, we again be- 
gin with a simple example of resonant photon storage in 
a homogeneously broadened atomic ensemble with neg- 
ligible spin- wave decay. It is important to note that, in 
contrast to the cavity model, the free space model gives 
access to many different spin-wave modes, which makes 
the retrieval efficiency dependent on how storage is car- 
ried out (see paper II). Therefore, optimization of storage 
alone is not a priori very practical. However, as shown 
in paper II, the optimization of storage alone is indeed 
useful because, in many cases, it also optimizes storage 
followed by backward retrieval. 

In order to have slightly simpler mathematical expres- 
sions, we work in the co- moving frame (see paper II), 
although the same argument can be carried out using 
the original time variable, as well. The complex number 
equations of motion on the interval t € [0, T] are then 
(see Ref. [l[ and paper II) 



diP{z, i) = -P(5, i) + iVd£{S, i) ^ 
dtS{S,i) = in{t)P{z,i), 

with initial and boundary conditions 

P(z~,0) = 0, 
5(z,0) = 0. 



(12) 

in{t)S{-z,i), (13) 
(14) 

(15) 
(16) 
(17) 



These equations are written using dimensionless vari- 
ables, in which (co-moving) time and Rabi frequency are 
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rescaled by 7 (t = <7 and 17 = 11/7) and the position 
is rescaled by the length L of the ensemble (z = z/L). 
£{z,t) describes the slowly varying electric field enve- 
lope, the input mode fin(i) satisfies the normalization 



constraint |£'in(f)| di = 1, d is the resonant optical 

depth, and VL{z) and Sinit) are for now assumed to be 
real. [To avoid carrying around extra factors of 2, d is de- 
fined as half of what is often referred as the optical depth: 
the steady-state solution with 17 = gives probe intensity 
attenuation \S{z = 1)|^ = e'^rf |£(z = 0)| .] The goal is 
to maximize the storage efficiency 



dz 



S{S,f) 



(18) 



with respect to 17 (t). A procedure analogous to that used 
in the cavity model in Sec. Ill Al yields equations of mo- 
tion (also referred to as the adjoint equations) for the 
Lagrange multipliers E{z,t), P(z,t), and S{z,t): 



di£ = iVdP, 

dlP = P + iVd£ + mS, 

diS = mP, 

with initial and boundary conditions 

£il,i) = 0, 
F(f,f) = 0, 
5(z,f) = 5(z,f). 



(19) 
(20) 
(21) 



(22) 
(23) 
(24) 



As in the cavity discussion in Sec. Ill Al these equations 
describe backward retrieval and provide a simple mathe- 
matical connection between optimality and time-reversal. 
In order to move up the gradient, one should update fl(t) 
according to 



1 



dz Im 



-P{z,i)S*(S,i) 



S*{S,i)P{z,i) 



(25) 



We showed in Ref. [l| and in paper II that, in the adi- 
abatic limit {Tdj ^ 1) and for a certain class of input 
modes of duration T ~ 1/(^7), one can achieve a univer- 
sally optimal (for a fixed d) storage efficiency that cannot 
be exceeded even if one chooses a different input mode. 
We showed that in that case the obtained control field 
will also maximize the total efficiency of storage followed 
by backward retrieval. However, this would not neces- 
sarily be the case for a general input mode in the non- 
adiabatic limit {Td'y ^1), which is precisely the limit, in 
which gradient ascent optimization becomes most useful. 
Moreover, for the case of forward retrieval, the control 
field that maximizes the storage efficiency does not max- 
imize the total efficiency of storage followed by retrieval 
even in the adiabatic limit. Thus, in Appendix [Cl we de- 
scribe how to use gradient ascent to maximize (still with 
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FIG. 4: Adiabatic (dotted) and optimal (solid) control fields 
for the storage followed by backward retrieval of a Gaussian- 
like input mode fin(t) (dashed) in the free space model with 
d = 10 and T = 5O/7 (a) and T = O.5/7 (b). Four optimal 
control pulses were obtained using four different initial guesses 
for the gradient ascent procedure. The adiabatic control field 
agrees with the optimal one in the adiabatic limit {Td'y 2> 1) 
(a) and deviates from it otherwise (b). 



respect to the storage control field) the total efficiency of 
storage followed by retrieval. 

As in the cavity model in Sec. Ill Al we now compare 
adiabatic shaping of the storage control field (see Ref. [l| 
and paper II) to the optimal shaping via gradient ascent. 
To compare with the results of paper II, we maximize the 
total efficiency of storage followed by backward retrieval 
rather than the storage efficiency alone. We assume that 
d — 10 and that £in{t) is the Gaussian-like input mode 
in Eq. (fTU]) . shown as a dashed line in Figs. H^a) and 
IDJb). We first consider the case T = 5O/7 and shape the 
storage control using adiabatic shaping (Sec. VI B of pa- 
per II) . Then we numerically compute the total efficiency 
of storage followed by complete backward retrieval using 
this storage control field (the total efficiency is indepen- 
dent of the retrieval control field provided no excitations 
are left in the atoms). The adiabatic storage control is 
shown as a dotted line in Fig. [H^a). Since for this input 
mode the adiabatic limit is satisfied {Td-f = 500 3> 1), 
the adiabatic storage control yields an efficiency of 0.66, 
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ih 100 




FIG. 5: The total efficiency of storage followed by backward 
retrieval for the Gaussian-like input mode in Eq. (|10p using 
adiabatic equations (dotted) and gradient ascent (solid) to 
shape the storage control field. The results are shown for the 
indicated values of d (= 1, 10, 100), as a function of Td'y. 
The dashed lines represent the maximum efficiency possible 
at the given d 



which is the maximum efficiency possible at this d 
For the same reason, the adiabatic control agrees with 
the control field computed via gradient ascent (solid line), 
which also yields an efficiency of 0.66. Fig. SJa) shows 
four solid lines (optimal control fields) corresponding to 
four initial guesses Vl{t)/^ = 0.2, 0.5, 1, and 1.5. As in 
the cavity model discussion in Sec. Ill A[ the difference be- 
tween the four optimal controls and the adiabatic control 
is inconsequential. 

Repeating the calculation for T = O.5/7, we obtain 
Fig. Hfb). Since the adiabatic limit {Td'y ^ 1) is no 
longer satisfied (^^7 = 5), the adiabatic approximation 
does not work and the adiabatic control differs from the 
optimal control and gives a lower efficiency: 0.24 vs. 0.58. 
As in Fig.lUJa), the four optimal control fields plotted cor- 
respond to different initial guesses il{t)/j = 1, 3, 5, and 
7. As in the cavity discussion, Fig. IDJb) indicates that, 
in the regime Td'y ~ 1, where the adiabatic approxima- 
tion no longer holds, the optimal control field acquires 
characteristics of the control field used in fast storage. 

As in the analysis of the cavity model in Sec. Ill A[ we 
now analyze the dependence on d and Td'y of the effi- 
ciency yielded by the adiabatic control shaping and the 
optimal control shaping. In Fig. [5l we compare the effi- 
ciency of storage followed by complete backward retrieval 
of the input mode in Eq. (fTU]) obtained using the control 
field shaped using the adiabatic equations (dotted lines) 
and using gradient ascent (solid lines). The efficiencies 
are plotted as a function of Td-y for three indicated val- 
ues of d {— 1, 10, 100). Horizontal dashed lines represent 
the maximum efficiency possible at the given d m. The 
dotted lines are the same as in Fig. 6(a) of paper 11. Sim- 
ilar to the corresponding discussion of the cavity model 
in Sec. Ill Al Fig. [5] confirms the predictions of papers II 
and III that efficient photon storage is not possible for 



Td'y < 1. It also illustrates that optimal control fields 
open up the possibility of cfRcient storage of input modes 
with a bandwidth that is an order of magnitude larger 
than the bandwidth allowed by the adiabatic storage. In 
addition, the same reasoning as in the cavity discussion 
leads to the conclusion that for this problem, gradient 
ascent most likely yields the global, rather than local, 
maximum at all values of Td'y. 

Various generalizations of the presented procedure can 
be made. First, the generalization to limited control 
pulse energy, (possibly reversible [s^l) inhomogeneous 
broadening, complex Q and 8^, and nonzero A, 7^, and 
7c can be carried out exactly as in the cavity case (Ap- 
pendix [B]). Second, in the case of backward retrieval, if 
the two metastable states are nondegenerate and have a 
frequency difference ujsg, one should incorporate an ap- 
propriate position-dependent phase shift of the spin wave 
of the form exp(— 2jAfcz), where Afc — Lujsg/c (see Sec. 
Vlll of paper II). Finally, another extension can be made 
for the cases when the total efficiency depends on the re- 
trieval control field (e.g. if 7s and/or 7c are nonzero). 
In those cases, one can simultaneously optimize with re- 
spect to both the storage and the retrieval control fields. 
However, one may then need to put a limit on the en- 
ergy in the retrieval control pulse since, for the case of 
7s ^ 0, for example, the faster one retrieves, the higher 
is the efficiency, and the optimal retrieval control field 
may, in principle, end up having unlimited power (e.g. 
an infinitely short tt pulse). 



III. OPTIMIZATION WITH RESPECT TO THE 
INPUT FIELD 

Although it is usually easier to optimize with respect to 
the control field, optimization with respect to the input 
mode can also be carried out in certain systems. For both 
classical and quantum light, the mode shape can often 
be controlled by varying the parameters used during the 
generation of the mode. For example, if the photon mode 
is created by releasing some generated collective atomic 
excitation, one can, under certain assumptions, generate 
any desired mode shape 1] . For the case of classical light, 
one can also shape the input light pulse simply using an 
acousto-optical modulator. An important advantage of 
optimizing with respect to the input mode is that the 
iterations can be carried out experimentally In t^iis 

section, we consider the maximization of light storage 
efficiency with respect to the shape of the input mode. 

The gradient ascent method, used in Sec.|TT]to optimize 
with respect to the control field, can be easily applied to 
the optimization with respect to the input mode shape 
both in the cavity model and in the free space model. 
Since one is interested in finding the optimal input mode 
shape, the optimization has to be carried out subject to 
the normalization condition ^"^ dt\£{a{t)\^ — 1. This con- 
dition can be included by adding an extra term with a 
Lagrange multiplier to the functional J to be optimized. 
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The iterations are then done as follows: one first inte- 
grates the storage equations for a trial input mode; then 
integrates the adjoint equations corresponding to back- 
ward retrieval (as in Sees. Ill A1 and HTB]) : then updates 
the trial input mode by adding to it a small correction 
proportional to the output of backward retrieval [~iP{t) 
in the cavity model or £{0,i) in the free space model]; 
and finally renormalizes the new input mode to satisfy 
the normalization condition. 

An important feature that distinguishes the optimiza- 
tion with respect to the input mode from the opti- 
mization with respect to the control field is the possi- 
bility of making finite (not infinitesimal) steps. Stan- 
dard gradient-ascent improvement [such as via Eqs. ([9|) 
and ()25|) ] is, in principle, infinitesimal due to its re- 
liance on the small parameter 1/A. Several decades ago, 
Krotov introduced and developed an important power- 
ful and rapidly converging global improvement method 
[Ij EH, m, 27, 28] that is not characterized by a small 
parameter. Largely thanks to the presence of the nor- 
malization condition on the input mode, this method can 
be applied to derive non-infinitesimal quickly converging 
updates for the problem of optimization of light storage 
efficiency with respect to the input mode. For the cavity 
model of Sec. Ill A[ this update is given by 

£Ut) ^ -^P{t), (26) 

followed by a renormalization of £in(i), while for the free- 
space model of Sec. IIIBI the update is given by 

£in(f)^f(0,?), (27) 

followed by renormalization. These updates precisely 
correspond to the time-reversal-based iterations sug- 
gested in Ref. [l| and explained in more detail in papers 
I and II. In these iterations, optimization of light stor- 
age with respect to the input field is done by carrying 
out storage of a trial input mode followed by backward 
retrieval, and then using the normalized output of back- 
ward retrieval as the input mode in the next iteration. 
The beauty of this update procedure is the possibility of 
carrying it out experimentally. In fact, the extension of 
this procedure to the optimization of storage followed by 
forward retrieval, suggested in Ref. [l| and in paper II, 
has already been demonstrated experimentally 0] • 

In the language of gradient ascent, one can still think 
of Eqs. (pS)) and (|27p as steps along the gradient. These 
steps are, however, finite, not infinitesimal. This allows 
one to think of time-reversal-based optimization with re- 
spect to the input mode as simple intuitive walk up the 
gradient. As shown in paper I, the fact that only one col- 
lective atomic mode can be excited in the cavity model 
makes the iterations of Eq. (pS)) converge to the opti- 
mum in a single step. Using the terminology of gradient 
ascent, the optimization with respect to the input field 
in the cavity model can, surprisingly, be achieved with a 
single large step up the gradient. 

We note that the optimization procedure discussed in 
this section can be easily generalized to include inhomo- 



geneous broadening and (for the case of Doppler broad- 
ened gases) the presence of velocity changing collisions. 
One can show that, even with these features, the iter- 
ative optimization procedure still works in exactly the 
same way by updating the input mode with the output 
of time-reversed retrieval. 



IV. OPTIMIZATION WITH RESPECT TO THE 
INHOMOGENEOUS PROFILE 

Having discussed optimization with respect to the con- 
trol field and the input mode, we now turn to the opti- 
mization with respect to the shape of the inhomogeneous 
profile. This optimization is most relevant in the con- 
text of controlled reversible inhomogeneous broadening 
(CRIB) [13]. The main idea of CRIB is that by intro- 
ducing inhomogeneous broadening into a homogeneously 
broadened medium (via Stark or Zeeman shifts, for ex- 
ample) and by optimizing the shape and width of this 
inhomogeneous profile, one can better match the absorp- 
tion profile of the medium to the spectrum of the in- 
coming photon mode and, thus, increase the storage ef- 
ficiency j30| . At the same time, one can minimize the 
losses caused by dephasing of different frequency classes 
with respect to each other by using an echo-like process 
triggered by a reversal of the inhomogeneousprofile be- 
tween the processes of storage and retrieval [l^, H^j ■ We 
refer the reader to papers I and III for a full list of ref- 
erences, and to Ref. [37| and paper III for examples of 
recent theoretical studies. 



A. Cavity model 

As in Sec. [TTl we begin the discussion with the theoret- 
ically simpler cavity model. Although one can, of course, 
optimize with respect to the inhomogeneous profile in the 
problem of storage alone (i.e. not followed by retrieval), 
in the context of CRIB it is more relevant to consider 
the problem of storage followed by retrieval with the re- 
versed inhomogeneous profile [30|. Moreover, although 
the approach can be extended to nonzero single-photon 
detuning and arbitrary control fields, we suppose for sim- 
plicity that the input mode £in{t) is resonant and that the 
storage and retrieval control pulses are tt pulses. Follow- 
ing the convention of Ref. [l| and papers I-III, we refer to 
this use of 7r-pulse control fields as fast storage and fast 
retrieval. 

We leave most of the mathematical details of the prob- 
lem to Appendix |DJ Here we only note that we describe 
the inhomogeneous profile by a discrete number of fre- 
quency classes. The index j labels the frequency class 
with detuning Aj from the center of the line containing 
a fraction pj of atoms (^jPj = !)• In Appendix iDl we 
show how to carry out optimization with respect to pj 
and/ or Aj. 
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FIG. 6: Comparison of the efficiency for storage followed by 
retrieval in the cavity model with and without controlled 
reversible inhomogeneous broadening (CRIB). We consider 
storage of the Gaussian-like input mode of duration T [Eq. 
(|10|) ] in a cavity with C = 50. (a) The figure shows the 
efficiency of fast storage followed by fast retrieval with a ho- 
mogeneous line (dash-dotted), fast storage followed by fast 
retrieval with a reversible optimized inhomogeneous profile, 
i.e. CRIB (circles), optimal storage and retrieval with a ho- 
mogeneous line as in Fig. [3ja) (solid), and the asymptotic 
value /{C +1)^ (dashed), (b) The optimal inhomogeneous 
width Ai for CRIB. 



We now present the results of gradient ascent opti- 
mization with respect to the inhomogeneous profile for 
a particular example. We suppose that the input pulse 
is the Gaussian-like mode in Eq. (fTU|) and that C = 50. 
The total efficiency of storage followed by retrieval, as a 
function of TC7, is shown in Fig. EUa) for various stor- 
age protocols. The dash-dotted line gives the efficiency 
of fast storage (i.e. storage obtained by applying a con- 
trol TT pulse on the \g) — |e) transition at the end of the 
input mode at time T) followed by fast retrieval using 
a homogeneous line. As discussed in papers 1 and 11, a 
homogeneous ensemble enclosed in a cavity has only one 
accessible spin-wave mode and can, therefore, fast-store 
only one input mode, which has duration T l/{C"f). 
As a result, the decay at TC'f » 1 of the efficiency repre- 
sented by the dash-dotted line is dominated by leakage of 
the input mode into the output mode and not by polar- 
ization decay. We now consider intoducing reversible in- 
homogeneous broadening and iteratively optimizing with 
respect to its shape (using Eq. (jD6p or Eq. (jPSp ). As ex- 



pected, the efficiency grows with each iteration indepen- 
dently of the choice of the number of frequency classes, 
the choice of A^, and the initial guess for pj. The land- 
scape in the control space, however, depends on the num- 
ber of frequency classes and on Aj . This landscape is also 
not as simple as in Sees. HIl and lllli i.e. there exist local 
maxima. We did not perform an exhaustive search, but 
out of all the initial configurations, number of frequency 
classes, and Aj distributions that we tried, the highest 
efficiencies were obtained for the cases when gradient as- 
cent converged to only two nonempty frequency classes 
with opposite detunings (we have not been able to come 
up with a simple physical reason for this). We therefore 
focus on the case of only two frequency classes with de- 
tunings ±Ai and optimize with respect to Aj [using Eq. 
(ID9p ]. The optimized efficiency is shown with circles in 
Fig. EJa). For TC7 less than about 0.75, it is optimal to 
have Ai = 0. For larger TC-f, the optimal Ai is shown 
in Fig. mb): at small TC-f, it scales approximately as 
oc {TC'y)^^ and then slower. The presence of two fre- 
quency classes and hence two accessible spin wave modes 
instead of one allows us to reduce the leakage error, so 
that the efficiency [circles in Fig. [H^a)] is now limited by 
polarization decay. 

Finally, we would like to compare the broadening- 
optimized efficiency to the homogeneous control- 
optimized efficiency. Repeating the optimization proce- 
dure of Sec. Ill AI for C = 50, we obtain the solid line in 
Fig. [6l[a) . The maximum efficiency possible at this C is 
C^/ (C-|- 1)^ and is shown as the dashed line. The dashed 
line and the solid line are the same as in Fig. [21 a) , ex- 
cept that now C — 50. The fact that the solid line in 
Fig. [6ja) lies above the circles indicates that we have not 
been able to identify any advantage of fast storage with 
CRIB compared to optimal storage in the homogeneous 
medium. Moreover, all inhomogeneous broadening con- 
figurations we tried to introduce into the optimized ho- 
mogeneous protocol converged back to the homogeneous 
profile. These results suggest that if one wants to store 
a single mode of known shape using a homogeneously 
broadened ensemble of A-type systems enclosed in a cav- 
ity and can shape and time the control field with sufficient 
precision, it may be better to use optimal homogeneous 
storage and not to use CRIB. 

It is, however, worth noting that we have only car- 
ried out the simplest optimization of fast storage with 
CRIB. In particular, the performance of fast storage with 
CRIB may be further enhanced by optimizing with re- 
spect to the time, at which the storage tt pulse is ap- 
plied. Such optimization represents an optimal control 
problem with a free terminal time Q and is beyond the 
scope of the present paper (although it can be carried 
out in a straightforward manner by repeating the opti- 
mization above systematically for different times of the 
TT-pulse application). 

It is also important to note that the use of CRIB in the 
cavity model may allow for implementing a multimode 
memory [l6j in the cavity setup. Unlike the free space 



11 



1 r 




FIG. 7: Comparison of optimized homogeneous-line storage 
with storage based on CRIB. For d = 100, the plot shows 
the efficiency of storage followed by backward retrieval of 
the Gaussian-like input mode of duration T [Eq. (|10p ] . The 
curves show results for fast storage and retrieval with a ho- 
mogeneous line (dash-dotted), fast storage and retrieval with 
an optimized reversible Gaussian (G) or Lorentzian (L) in- 
homogeneous profile, i.e. CRIB (solid lines labeled G and L), 
storage and retrieval with a homogeneous line using adiabatic 
(dotted) or optimal (unlabeled solid line) control field shaping 
(same as in Fig. [5]l, and the asymptotic value (dashed). 



model, which allows for the storage of multiple temporal 
input modes using, e.g., Raman- or EIT-based protocols 
[JH, m, nil , the homogeneously broadened cavity model 
only has a single accessible spin-wave mode. Therefore, 
if we do not use CRIB or some other inhomogeneous 
broadening mechanism, it can only store a single input 
mode 



B. Free space model 

Having discussed the optimization with respect to the 
inhomogeneous profile in the cavity model, we note that 
the same procedure can be carried out for the free space 
model in an analogous manner. The appropriate update 
equations are listed at the end of Appendix [Pl 

In paper III, we compared storage using adiabatic con- 
trol shaping in a homogeneous ensemble to fast storage 
with CRIB. We found that fast storage with CRIB can 
indeed do better than adiabatic homogeneous storage for 
Td^ ^ 1. We show now that this result was mainly 
due to imperfect control field optimization outside of the 
adiabatic limit and that, in the present work, we have 
not been able to identify any advantages of fast storage 
with CRIB compared to optimal homogeneous storage. 
We consider storage of the resonant Gaussian-like input 
mode in Eq. (|10p in a free space atomic ensemble with 
d — 100 followed by backward retrieval. The total effi- 
ciency for various storage protocols is shown in Fig. [7] as a 
function of Td'y. The dash-dotted line and the two solid 
lines labeled G and L are taken from Fig. 8(a) of paper 
III. The dash-dotted line is the efficiency of fast storage 



followed by fast backward retrieval using a homogeneous 
line. The two solid lines labeled G and L are obtained 
using fast storage with optimal- width reversible inhomo- 
geneous broadening with Gaussian profile and Lorentzian 
profile, respectively. Although the optimization with re- 
spect to the inhomogeneous width can be done efficiently 
via gradient ascent [using Eq. (IDlip j. we have already 
performed this optimization in paper III by sampling a 
sufficiently large set of inhomogeneous widths. The re- 
maining third solid line and the dotted line (both taken 
from Fig. [5]) correspond to homogeneous storage with op- 
timal storage controls (solid) and with adiabatic controls 
(dotted). The dashed line (also from Fig. [5]) is the max- 
imum possible efficiency at this d. The plot shows that 
while adiabatic control field shaping (dotted) makes ho- 
mogeneous storage less efficient for some values of Td'y 
than fast storage with CRIB (solid lines labeled G and L), 
optimal control field shaping (unlabeled solid line) may 
enable homogeneous storage to be more efficient than fast 
storage with CRIB at all values of Td^. 

As in Sec. IIV A| we note, however, that we have pre- 
sented only the simplest optimization of CRIB and that 
the full investigation of the advantages of CRIB is be- 
yond the scope of the present paper. In particular, the 
CRIB efficiency may be enhanced by optimizing with re- 
spect to the time, at which the storage 7r-pulse is applied. 
Moreover, CRIB might be useful in circumstances such 
as when a homogeneously broadened three-level system 
is not available, when more complicated inputs (such as 
time-bin qubits) are used, or when precise shaping and 
timing of the control pulse is harder to achieve than con- 
trolled reversible broadening. Finally, CRIB-based mem- 
ories may even be implemented without any optical con- 
trol fields [13]. 



V. SUMMARY 

In conclusion, we have shown that the powerful numer- 
ical optimal control method of gradient ascent allows one 
to obtain simple intuitive understanding and to achieve a 
significantly improved efficiency and a higher bandwidth 
in the problem of photon storage in A-type atomic en- 
sembles. First, we showed how to apply gradient ascent 
to numerically compute optimal control fields even out- 
side of the adiabatic limit both with and without a con- 
straint on the energy in the control pulse. In particular, 
this opens up the possibility of efficient storage of in- 
put modes that are an order of magnitude shorter (and 
hence an order of magnitude larger in bandwidth) than 
the shortest modes that can be efficiently stored using 
adiabatic control field shaping. Second, we showed that 
gradient ascent provides an alternative justification for 
the often discussed connection between optimality and 
time-reversal in photon storage, as well as for the itera- 
tive time-reversal-based optimization procedure with re- 
spect to the input field suggested in Ref . [J] , discussed in 
detail in papers I, II, and III, and demonstrated experi- 
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mentally in Ref. In particular, we confirmed that the 
iterative procedure works even in the presence of inhomo- 
geneous broadening and (for the case of Doppler broad- 
ened gases) in the presence of velocity changing collisions. 
Finally, we showed how to use gradient ascent to optimize 
with respect to inhomogeneous broadening and demon- 
strated how this can significantly increase the efficiency 
of fast storage followed by fast backward retrieval in the 
presence of controlled reversible inhomogeneous broad- 
ening (CRIB) [s^l- Provided one is interested in storing 
a single input photon mode of known shape and provided 
the control pulses can be generated with sufficient pre- 
cision, we have not, however, been able to identify any 
advantages of CRIB-based photon storage compared to 
photon storage with optimal control pulses in homoge- 
neously broadened media. 

In general, gradient ascent methods do not guarantee 
the attainment of the global maxima. The global maxi- 
mum is, however, indeed attained for our problem in the 
regimes where this maximum is known. This strongly 
suggests that, for the optimization with respect to the 
input mode and with respect to the storage control, gra- 
dient ascent may indeed be yielding the global optimum. 
We also note that one can optimize simultaneously with 
respect to various combinations of the control param- 
eters simply by simultaneously updating each of them 
along the corresponding gradient. One can also include 
other possible control parameters that are available in 
a given experimental setup but have not been discussed 
in the present paper. For example, for the case of pho- 
ton storage in solid-state systems, one can consider opti- 
mizing with respect to the number of atoms put back 
into the antihole [s^, [s^ or with respect to a time- 
dependent reversible inhomogeneous profile. Other light 
storage systems, such as photonic crystals [i^ or cavity 
models where the cavity field cannot be eliminated, are 
also susceptible to gradient ascent optimization. There- 
fore, we expect the optimization procedures described in 
the present paper to allow for increased efficiencies and 
increased bandwidths in many current experiments on 
quantum memories for light, many of which are narrow- 
band and suffer from low efficiencies. Such improvements 
would facilitate advances in fields such as quantum com- 
munication and quantum computation. 
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APPENDIX A: DERIVATION OF THE ADJOINT 
EQUATIONS OF MOTION IN THE CAVITY 
MODEL 

In Sec. Ill Al we omitted the derivations of the adjoint 
equations of motion ^ and ([5]) and of the correspond- 
ing boundary conditions ^ and ([7]). We provide these 
derivations in this appendix. 

Varying J given in Eq. ([3]) with respect to S, S* , P, 
and P*, we obtain 

6J = S{T)6S*{T) 

+ J dtS* (-SS + iflSP^ 

+ c.c, (Al) 

where "c.c." means complex conjugate taken of the whole 
expression after the equal sign. Integrating by parts the 
terms containing time derivatives, we obtain 

dJ = S{T)6S*{T) 

-~P*{T)SP{T)+ [ dtP*SP+ [ dtp*(inss 
Jo Jo ^ 

-7(1 + C)SP) 

~S*{T)SSiT)+ [ dtS*SS+ [ dtS*{inSP) 

Jo Jo 
+ c.c. (A2) 

Since the initial conditions are fixed, we have here used 
6S{0) = <5P(0) = 6S*{0) = SP*{0) = to simplify the 
expression. 

The optimum requires that SJ — for any variations 
6P and SS. Hence we collect the terms multiplying, e.g., 
SP{T) and set the result to zero. Carrying out this proce- 
dure for 6P{T), SS{T), and their conjugates, we obtain 
the boundary conditions ^ and ([7]). Collecting terms 
proportional to SP, SS, and their conjugates, we obtain 
adjoint equations of motion ^ and ([5]). 

APPENDIX B: CONTROL FIELD 
OPTIMIZATION IN THE CAVITY MODEL: 
GENERALIZATION 
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In Sec. Ill Al we showed how to perform control field 
optimization in the simplest possible version of the cav- 
ity model: a resonant input mode with a real envelope 
was stored using a control pulse with a real envelope and 
unlimited power into a homogeneously broadened ensem- 
ble with infinite spin-wave lifetime. In this appendix, 
we show how to optimize the control field in a more 
general model that includes the possibility of complex 
control field envelopes fl(t) and input mode envelopes 
Sinit), nonzero single-photon detuning A and spin wave 
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decay rate 7s, and (possibly reversible [30|) inhomoge- 
neous broadening such as Doppler broadening in gases 
or the broadening of optical transitions in solid state im- 
purities [sg'I . For the case of Doppler broadened gases, we 
also include velocity changing collisions with rate 7c. We 
also show how to take into account possible experimental 
restrictions on the strength of the classical control fields. 

Using the notation of paper III, the complex number 
equations describing the generalized model are 

P, = -[7 + z(A + A,)]P,~7C^F + zf^5, 

+zv/2^^fin + Ici^jP - Pj), (Bl) 
Sj = --f,S,+in*P,+-/ci^S^Sj), (B2) 

where j labels the frequency class with detuning Aj from 
the center of the line containing a fraction pj of atoms 
C^jPj = 1) and where the total optical and spin polar- 
izations are P = J2k \fPkPk and S = ^/pkSk, respec- 
tively. The terms proportional to 7c describe completely 
rethermalizing collisions with rate 7c [4l|. One can, of 
course, also take 7c to be different for P and S. For ex- 
ample, if 7c 7, which is often the case, one can drop 
the terms proportional to 7c in Eq. (|B1|) [i^ . In addi- 
tion to moving atoms from one frequency class to the 
other, collisions also result in line broadening, which can 
be taken into account by increasing 7 [i^]. We assume 
that the goal is to maximize the efficiency 77s = [^(T)^ 
of storage into the symmetric mode S{T) with respect to 
the control pulse n{t) for a given input mode shape £i-a{t) 
satisfying the normalization condition jj" dt \£\i^(t)^ = 1. 
A procedure very similar to that described in Sec. Ill Al 
and in Appendix [X] yields the following equations of mo- 
tion for the adjoint variables: 

Pj = [^-i{A^Aj)]Pj+^C^jP^inSj 

-1c{Vp]p-p,), (B3) 



isSj + in*Pj - -fc{^/PjS - Sj), 



(B4) 



where P ^ J2k y/P^Pk and S ^ J2k \/PkSk- The cor- 
responding initial conditions for backward propagation 
are 



PAT) 
SAT) 



0, 



p-S{T). 



(B5) 
(B6) 



Similarly to Sec. Ill Al after taking an initial guess for 
n{t) and solving for Pj, Sj, Pj, and Sj, one updates n{t) 
by moving up the gradient 



nit)^n{t) + jzJ2{s*Pj-PjS* 



(B7) 



Short input modes and/or large single-photon detun- 
ing A require control pulses with large intensities that 
might not be available in the laboratory. There exist 
ways to include a bound on the control field amplitude 
0. Alternatively, one may want to consider a slightly 



simpler optimization problem with a limit on the control 
pulse energy \fl(t)\^ dt < E for some E [43|. In order 
to carry out the optimization subject to this constraint, 
one should first carry out the optimization without the 
constraint and see whether the optimal control satisfies 
the constraint or not. If it does not satisfy the constraint, 
one has to add a term ^'{E — 1^^(01^ dt) to J, so that 
the update becomes 



nit)^n{t) + \ 

A 



iJ2{s*p,-p,s*)^fi'n{t) 



(B8) 



where fj,' is adjusted to satisfy the constraint. By redefin- 
ing fi' and A, this update can be simplified back to Eq. 
(|B7|) followed by a renormalization to satisfy the con- 
straint. Depending on how severe the energy constraint 
is, one can then sometimes (but not always) further sim- 
plify the update by completely replacing n{t) with the 
gradient [i.e. set X = fi' in Eq. (|B8p ] followed by a renor- 
malization of fl{t), as is done, for example, in Ref. [l3| 
for the problem of laser control of chemical reactions. 

We note that these optimization protocols can be triv- 
ially extended to the full process of storage followed by 
retrieval, which, in the presence of inhomogeneous broad- 
ening, one might not be able to optimize by optimizing 
storage and retrieval separately. Similarly, one may in- 
clude the possibility of reversing the inhomogeneous pro- 
file between the processes of storage and retrieval [30|. 



APPENDIX C: CONTROL FIELD 
OPTIMIZATION IN THE FREE-SPACE MODEL: 
GENERALIZATION TO STORAGE FOLLOWED 
BY RETRIEVAL 



In Sec. IIIBI we showed how to use gradient accent 
to find the control field that maximizes the storage effi- 
ciency. However, the obtained storage control field does 
not always maximize the total efficiency of storage fol- 
lowed by retrieval. Therefore, in this appendix, we con- 
sider the maximization of the total efficiency of storage 
followed by retrieval with respect to the storage control 
field. While we demonstrate the procedure only for the 
case of forward retrieval, the treatment of backward re- 
trieval is analogous. 

We suppose that the control field flit) consists of a 
storage control pulse on t G [0, T] and a retrieval control 
pulse on t G [Tr, T{]. We want to optimize with respect to 
the former given the latter and the input mode (note that 
the total efficiency is independent of the retrieval control 
for sufficiently strong retrieval control pulses, and it is 
therefore often less important to optimize with respect 
to the retrieval control pulse). Here < T < T,- < T{, 
and the subscripts in and Tf stand for "retrieval" and 
" final" . The time interval [T, T^] corresponds to the wait- 
ing (i.e. storage) time between the processes of storage 
(which ends at t = T) and retrieval (which begins at 

t = Tr). 
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We suppose that storage is described by Eqs. (|T2|l - (fT7|) 
on t G [0,T]. Then forward retrieval that follows after 
the storage time interval [T, T^] is described by the same 
equations ((T2])-(IT1]) but on the time interval t e [Tr,Tf] 
with initial and boundary conditions 



£{0,i) = 0, 

P{S,fr) = 0, 



(CI) 
(C2) 
(C3) 



where = T^-f (similarly, ff = Ta)- Eq. assumes 
that the polarization has sufficient time to decay before 
retrieval starts, while Eq. (IC3P assumes that spin- wave 
decay is negligible during the storage time. The goal is 
to maximize the total efficiency of storage followed by 
retrieval. 



Vtot 



di\£{l,i)f . 



(C4) 



with respect to the storage control field. Constructing J 
and taking appropriate variations, we obtain initial and 
boundary conditions for backward propagation: 



£{l,i) = £(l,f) for i e [T;,Tf], (C5) 

P(5,ff) = 0, (C6) 

S{S,ff) = 0, (C7) 

= for ie [0,t], (C8) 

P(i,f) = 0, (C9) 

S{S,f) ^ S{S,f,). (CIO) 



and 



By taking the variational derivative of J with respect to 
f2(t) on the storage interval, we find that the update is 
exactly the same as for the optimization of storage alone 
and can be done via Eq. (^5)1 . 

We note that if the retrieval control pulse leaves no 
atomic excitations, one can obtain the same optimization 
equations by solving the storage optimization problem in 
Sec. lIIBl but changing the function to be maximized from 
the number of spin-wave excitations dzS{z, T)S* (z, T) 
to the complete retrieval efficiency from S{z, T) [Eq. (6) 
of Ref. [l|]. It is also worth noting that the derivation 
presented here can trivially be extended to apply to back- 
ward (instead of forward) retrieval and to include com- 
plex Q, and £in, (possibly reversible [s^]) inhomogcncous 
broadening, and nonzero A, 7s, and 7c. 



APPENDIX D: OPTIMIZATION WITH RESPECT 
TO THE INHOMOGENEOUS PROFILE: 
MATHEMATICAL DETAILS 

In Sec. IIVI we presented the results on the optimiza- 
tion of photon storage with respect to the inhomogeneous 



broadening without providing the mathematical details. 
In this appendix, we present these details. 

We first consider the cavity model, but turn briefly to 
the free-space model at the end of this appendix. We 
suppose for simplicity that the input mode £in{i) is res- 
onant and that the storage and retrieval control pulses 
are tt pulses at t = T and t = T-^, respectively. In order 
to simplify notation, we define Xj = ^Jp], satisfying the 
normalization x'^ ~ 1. The storage equation (|B1[) on 
the interval t E [0, T] then becomes 

Pj = -{^ + iAj)P,-jCxjP + i^/2^Xj£in,im 

with P = XkPk and with the initial condition Pj{0) = 
0. A TT-pulse at t — T mapping P onto S followed by an- 
other TT-pulse at t = Tr mapping S back onto P result in 
an overall 2tt pulse, so that PjiT^) = —Pj{T). Assuming 
the broadening is reversed at some time between T and 
Tr, the equations for retrieval on the interval t G [Tr,Tf] 
are 

Pj = -{-f -iAj)Pj -"fCxjP. (D2) 
The total efficiency of storage followed by retrieval is then 



?7tot 



Tt 



dt\£outit)\' 



Tt 



dt 



iy/2^'P{t) . (D3) 



One can show that the equations of motion for the 
adjoint variables (i.e. the Lagrange multipliers) Pj are 



P, = h + iA,)P, 



for t e [Tr,Tf\ with Pj{Tf) 

P, = ij-iAj)Pj+jCx,P 



--fCxjP-2jCxjP (D4) 
and 



(D5) 



for t e [0,T]_ with Pj{T) = -Pj{Tr), where we defined 
P — J2k^kPk- The last term in Eq. (|D4p describes an 
incoming field that is the time-reverse of the retrieved 
field. Assuming we are optimizing with respect to Xj, 
the update is 



(D6) 



followed by a rescaling of all Xj by a common factor to 



ensure the normalization x'j 



A. 



1 . Here Aj is given by 







7CRe 


[ dt+ [ dt 




Jo JTr 



{p*p 



p*p,) 



-V27CIm / dt£i,,P* 
Jo 

+2-fCRe J dtP*P, 



(D7) 



where Re denotes the real part. Numerics show that the 
update can usually be simplified in a way that avoids 
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the search for convenient values of A and does not lose 
convergence. Specifically, taking A in Eq. (|D6p . we 
obtain 



(D8) 



followed by renormalization. By defining a particular 
functional form for the dependence of Xj on Aj , one could 
also consider optimization with respect to only a few pa- 
rameters, such as, for example, the width Aj and the 
degree of localization n of the inhomogeneous profile of 
the form pj = cx 1/ [1 + (Aj/Ai)"]. 

Equivalently, instead of optimizing with respect to Xj , 
one can optimize with respect to Aj. To illustrate this 
procedure, we consider a simple optimization procedure 
with respect to a single parameter, the inhomogeneous 



By adjusting fj and a;^, one can choose a particular in- 
homogeneous profile shape (e.g. Lorentzian, Gaussian, or 
a square) and optimize with respect to its width. 

Having discussed the cavity case, we now list the cor- 
responding free-space results. In free space, the update 
of Xj via Eq. (jD8|l would use 









/ dz 


[ di + 


/ dt 


'o 


Jo 





{P*£ + £*Pj) . 

(DIO) 



width Ai. We write A^- = Aifj for some fixed dimension- Similarly, the update of Aj = A1/7 would be 
less parameters fj and consider maximizing the efficiency 
with respect to Aj for fixed Xj and fj. The equations of 
motion and the initial conditions for both Pj and Pj stay 
the same as in the optimization with respect to Xj while 

the update becomes - , 1 , 



At 



dt 







f dz 


/ di- di 


JO 


Jo Jf, 



(Dll) 
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